Why is lead so kinky? 
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We revisit the problem of the kink in the charge radius shift of neutron-rich even isotopes near the 
TV — 126 shell closure. We show that the ability of a Skyrme force to reproduce the isotope shift is 
determined by the occupation of the neutron lin/2 orbital beyond TV = 126 and the corresponding 
change it causes to deeply-bound protons orbitals with a principal quantum number of 1. Given 
the observed position of the single-particle energies, one must either ensure occupation is allowed 
through correlations, or not demand that the single-particle energies agree with experimental values 
at the mean-field level. 
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PACS numbers: 21. 10. Ft, 21.60.Jz, 21.30.-x, 27.80.+W 

The evolution of charge radii across the isotope chart 
is one of the most basic nuclear structure observables, as 
it provides a particularly useful characterisation of the 
proton distribution that can be accessed by a variety of 
experiments [l[ . We define the charge radius isotope shift 
as the differences between the mean squared charge ra- 
dius, (r^), of a series of isotopes and that of a given ref- 
erence isotope ( 208 Pb for lead and 210 Po for polonium). 
Theoretically, charge radii have been traditionally stud- 
ied within the droplet model Q, which captures quali- 
tatively their evolution across the nuclear chart. Yet, in 
some specific cases, quantum shell effects dominate the 
density distribution and provide a departure from smooth 
systematic trends. Perhaps the most well-known exam- 
ple of these abrupt changes is the kink in the isotope shift 
of even lead isotopes as one passes through the TV = 126 
magic number. 

A summary of the experimental results [H, H| around 
this neutron shell closure is given in Figures [1] and [2] for 
lead and polonium, respectively. For lead, the fitted solid 
line is of the form: 
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Linear regression gives 



mi = 0.0598(7) fm 2 
m 2 = 0.1203(9) fm 2 , 

i.e. the slope of the shift is observed to double at 
A = 208 to a very good approximation. This abrupt 
change in change radius cannot be explained within the 
droplet model [2|. Recent experimental work using laser 
spectroscopy techniques have identified a similar kink in 
neutron-rich polonium isotopes above the Z = 82 shell 
closure 3 (see also Fig. [2} in agreement with older radon 
and radium data 0]. Moreover, the details of the changes 
of proton and neutron radii in isotopes around 208 Pb have 
also gained renewed interested in the form of the neutron 
skin, which correlates strongly with nuclear matter prop- 
erties [H, 0]. On the neutron deficient side, the onset of 
deformation is also probed by measurements of isotope 



shifts 0]. In this letter, we propose a new mechanism to 
explain the existence of the kink in lead isotope shifts us- 
ing density functional calculations supplemented by pair- 
ing effects. 

On the theoretical side, mean-field models, or equiv- 
alently density functional theories have been widely ap- 
plied to the systematic study of all observed and hypoth- 
esized nuclear isotopes Q. Both the Skyrme-Hartree- 
Fock (SHF) and the Relativistic Mean-Field (RMF) ap- 
proaches are able to give broadly good descriptions of 
many nuclear properties across the nuclear chart, includ- 
ing radii |{|. In particular, shell effects are naturally in- 
cluded in such quantum mechanical calculations. The 
kink in the isotope shift for lead, however, is still a some- 
what challenging observable, for which a full theoretical 
understanding remains elusive Note, in particular, 
that beyond mean-field correlations have a very small ef- 
fect on the charge radii in this area, close to the shell 
closure This suggests that any missing theoreti- 

cal component should have a rather general origin, most 
likely independent of the particular microscopic picture. 

The ability of different mean-field methods to repro- 
duce (or not) the isotope shift in lead has often been used 
to benchmark different models 11 - 13| . Original Skyrme 
parameterizations were unable to reproduce the isotope 
shift in lead, whereas RMF parameterizations seemed to 
be able to do so. Since the early 1990s, the situation 



has been studied by several authors I10h121 Il4j . Sharma 
et al. [lH and Reinhard and Flocard [l2| took the hint 
that the spin-orbit mean-field resulting from the SHF and 
RMF approaches has a different isotopic dependence [lj| . 
Reinhard and Flocard extended the Skyrme approach 
with an extra degree of freedom in the spin-orbit channel 
[l2| . In their notation, the Energy Density Functional 
(EDF) due to the spin-orbit interaction can be written: 
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with p denoting the total particle density, p q , the proton 
(or neutron) density, J, the spin-orbit current for both 
types of particle and J q , that for a particular type. 
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FIG. 1. Isotope shifts are given by the difference in the mean 
square charge radius between a series of even isotopes, de- 
noted by their mass number A, and that of 208 Pb. Across 
the N = 126 shell closure, a strong increase in the slope of 
the experimental data (diamonds) is observed. Theoretical 
predictions, obtained with different Skyrme parametrizations, 
are also presented. Only a handful of these Skyrme sets are 
able to reproduce the increase of slope above N = 126. 



With the choice 64 = 6 4 , one retrieves the form from 
the original Skyrme force as posited by Bell and Skyrme 



171 Il8j , while for b' 4 = one has an RMF-type spin-orbit 



mean field. One need not apply either of these rules, but 
rather allow a free fit of both 64 and b' 4 . Indeed, the best 



121 ] - the widely-applied force SkI4 - had 



X 2 found in Rcf. 
64 ~ — 6 4 . Furthermore, it was found that a qualitative 
description of the kink could be found with the original 
Bcll-Skyrme spin-orbit functional if one included the iso- 
topic shift data in the fit. This, however, can only be 
achieved at the expense of the overall fit quality through 
an unusually low effective mass, as is the case for force 
SkI5 [13 • 

The usual explanation given for the reproduction (or 
lack of) of the kink is the position of the 2g 9 / 2 neutron 
orbital above the N=126 shell gap 12, Ijj], which in turn 
is particularly sensitive to the spin-orbit component of 
the EDF. If this orbital is subject to a weak spin-orbit 
force, it will not be so deeply-bound and its radius will 
be commensuratcly larger. Consequently, the pull on the 
proton states through the symmetry energy will be some- 
what pronounced and might be able to explain the sud- 
den change in charge radius. It has also been noted that 
the l«ii/2 neutron state is sometimes occupied through 
pairing (typically BCS) probabilities. Yet, the Hn/ 2 or- 
bital is consistently ~1 fin smaller in radius than the 
2.99/2 [I2]s so a description based on the radius of oc- 
cupied neutron orbitals cannot give a full explanation. 
In this letter, we seek to clarify the important role of 
the occupation of the neutron lin/2 orbital in producing 
the correct isotope shift. We advocate that the overlap 
between neutron and proton orbitals with the same prin- 
cipal quantum number, n = 1, plays a major role. 



FIG. 2. The same as Fig. [T] for polonium. The reference 
isotope is 210 Po. 



The motivation for this study comes partly from a re- 
cent evaluation of most known Skyrme parameterizations 
for their ability to fit nuclear matter properties [2(J . We 
would like to see how those forces that passed a series of 



nuclear matter constraints perform in finite nuclei [21 



Note that the spin-orbit term of the mean field is inac- 
tive in infinite nuclear matter. Consequently, functionals 
which perform well in infinite systems might not be able 
to reproduce observables like the isotope shift. A slight 
subset of the "good" forces is used in the present work 
to serve as an example. They are: SKRA [22j, fitted to 
a realistic nuclear matter equation of state (EoS) along 
with some finite nuclear properties; and NRAPR [23[, fit- 
ted to the EoS, with adjustment of the spin-orbit force to 
optimise binding energy and radii in some doubly-magic 
nuclei. We also consider the well known forces SLy4 [24| , 
a widely-used parametrization from the Lyon group; and 
the already mentioned SkI4 12j, explicitly adjusted to 
reproduce isotope shifts in lead. We note that only one 
of these parameter sets, SkI4, uses the extended spin- 
orbit force of Eq. ([2]). In addition, we re-adjust SLy4, as 
stated below, to analyse the case of the extended form. 
We would also like to stress that our conclusions are gen- 
eral and do not depend on the specific Skyrme functional 
at hand. 

Let us briefly review the spin-orbit mean field prop- 
erties of some of the Skyrme models used here. As 
it stands, the spin-orbit parameter for NRAPR, Wo = 
2b A = 2b' A = 41.958 MeV fm 5 is very low (typical val- 
ues for other forces are Wo ~ 120 MeV fm 5 ). This 
results, in particular, in a proton shell-gap so small at 
Z=82 as to destroy the doubly-magic nature of the nu- 
cleus. We therefore replace NRAPR by constructing a 
modified force, NRAPRii, in which the value of the Wo is 
doubled, leaving other parameters unchanged. Further, 
to analyse the effect of a RMF-type spin-orbit force, we 
introduce a modified SLy4 force, with the original pa- 
rameters (64 = b' A = 61.5 MeV fm 5 ) changed to 64 = 75.0 
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FIG. 3. Neutron single-particle energies around the Fermi 
surface in 210 Pb for five sets of Skyrme forces. States with a 
significant BCS occupation, > 3%, above N=126 are embold- 
ened . Whenever the lin/2 state is substantially populated, 
the kink in isotope shift can be reproduced. 



MeV fm 5 and b' 4 = 0.001 MeV fm 5 . The change in 64 
for SLy4mod was designed to keep the magnitude of the 
spin-orbit field in 208 Pb the same as SLy4, while alter- 
ing only its isovector properties. This only modifies total 
energies and radii of closed shell nuclei within 2%. 

The calculated isotope shifts for these five Skyrme pa- 
ramctcrizations in lead and polonium arc shown in Figs.[T] 
and Fig. [5J respectively. Calculations were performed 
with a spherical Hartree-Fock code, including BCS pair- 
ing with a delta volume pairing force Q . Deformed calcu- 
lations for 186 Pb, not included here for brevity, indicate 
that all these forces also provide a plausible description 
of non-spherical isotopes. SLy4mod, NRAPRii and SkI4 
are able to produce a sizable kink in the isotope shift. 
Note, however, that SkI4 does so by construction. Most 
Skyrme forces, including SLy4 and SKRA, struggle to re- 
produce the kink. The forces which have a kink in lead 
do produce one in polonium. The mechanism underly- 
ing the kink should therefore be quite general within this 
mass region. 

A key difference between the forces showing the kink 
and the others is illustrated in Figure This shows the 
neutron single-particle levels in 210 Pb, which are only 
slightly rearranged from those in 208 Pb. In particular, 
the holding of the levels above the N=126 shell closure 
indicates those orbits outside the 208 Pb core with a BCS 
occupation above 3% (25|. A key point of this Letter 
is the observation that the lixi/2 state is substantially 
occupied only in the forces which reproduce the kink. 
The occupation of this state is vital for getting a kink in 
other isotopic chains. 

The effect of the lin/2 occupation on the charge radii 
is explored in Figure |4j Each frame shows the radius 
isotope shift of every individual occupied proton orbitals 
from the lowest ls 1 / 2 state up to the l/in/2 orbital below 
the Z = 82 shell gap. For the forces in which the H11/2 
neutron orbital has a significant occupation above N = 



126, there is an increase in the isotope shift of the proton 
states with principal quantum number n= 1. States with 
higher principal quantum numbers show a more mixed 
behaviour, n = 2 states have less marked differences 
across forces, while the parameterization dependence of 
the 3si/2 state is concentrated on the A < 208 side. The 
combined contribution of the n = 1 states explains the 
presence of the kink in the isotope shift 14 1 . 

From a nuclear matter perspective, one can say that 
the strong nuclear symmetry energy acts to increase the 
overlap between all the proton states and the overall nu- 
clear density. The effect will be enhanced if the overlap 
between wave functions is maximal, which occurs when- 
ever the nodal structure is the same. The latter is deter- 
mined by the principal quantum number. To test this, we 
calculate the radial overlaps in 208 Pb between the (un- 
occupied) lin/2 an d 2g 9 /2 neutron orbitals with each of 
the occupied proton orbitals. Since the particles in each 
orbital interact via the nodeless mean field, and because 
of the symmetry energy tending to favour overlapping 
states, the radial overlaps serve as a useful proxy for un- 
derstanding the isotope shift. These arc shown in Figure 
[5] for NRAPRii, but similar results hold for the other 
Skyrme forces. A clear dominance of the H11/2 overlaps 
over the 2g 9 / 2 ones is seen. Hence, when neutrons are 
added to the Kn/2 state, the proton states are attracted 
to larger radii to maximally overlap with the extra neu- 
trons. 

This leaves the question of the observed locations of 
the two neutron orbitals in question. Experimentally, 
the 2gq /2 orbital is more deeply bound than the lin/2 or- 
bital [26(| , though separated only by 1 MeV. SLy4mod, 
SkI4 and NRAPRii predict a smaller (or even a nega- 
tive) splitting between these orbitals. If one identifies 
the mean-field single-particle energies with experimen- 
tally measured values, the mechanism we propose can 
only work if the pairing interaction allows for the scat- 
tering of enough Cooper pairs of neutrons into the lin/2 
orbital. In general, however, mean-field single-particles 
energies should not necessarily agree with measured nu- 
clear energy levels 27 1. 

As previously observed, forces with an RMF-stylc spin- 
orbit force, like SkI4, or a low effective mass, also tend 
to give a level ordering favouring the occupation of the 
1*11/2 neutron orbit. According to our explanation, this 
should lead to a good reproduction of the lead isotope 
shifts. Our observation also gives scope for a) forces in 
which energy levels themselves are fitted, thus leading 
to effective masses close to m*/m = 1 [26| and provid- 
ing correlations that allow for a sufficient occupation of 
the H11/2 orbital, or b) a djust ments of the single-particle 
levels via the tensor force [28| to reproduce correctly this 
benchmark, without damaging good nuclear matter prop- 
erties. 

In conclusion, we have demonstrated that the repro- 
duction of the isotope shift in lead is by and large deter- 
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FIG. 4. Lead isotope shifts of individual proton orbitals for five sets of Skyrme forces. Filled squares correspond to SLy4, open 
squares to SLy4mod, filled circles to SKRA, open triangles to SkI4 and filled triangles to NRAPRii. The results highlighted in 
red are those associated to forces which reproduce the isotope shift. Above the shell closure, the n — I states in the leftmost 
column contribute substantially to the kink. 
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that are deeply-bound. This provides a larger "pull" of 
neutron states on proton orbits via the symmetry energy 
and allows for the reproduction of the well-known kink 
across the N = 126 shell gap. Our explanation is in con- 
trast to the more traditional account, based on the effect 
of the 2<7g/ 2 on the proton density. For our calculations, 
we have tweaked the spin-orbit component of a few mean- 
field parametrizations to lower the lin/2 single-particle 
energy and increase its occupation. The physical origin 
of this occupation, however, is not important for the gen- 
eral mechanism presented here. 

This work has been supported by STFC grants 
ST/J00051/1, ST/J500768/1, and ST/1005528/1. 



FIG. 5. Radial overlaps between all the occupied proton or- 
bitals in 208 Pb and the neutron 2g 9 /2 (left bars) and lin/2 
(right bars) states computed with the NRAPRii parametriza- 
tion. An overall predominance of overlaps where both states 
have n — 1 is found. 



mined by the occupation of the lin/2 neutron orbital. 
Since this is an n = 1 orbital, it overlaps more strongly 
with the majority of the proton orbitals, including those 
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